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SU»MARY 

A nonlinear analysis which includes the effect of^mismatch and filler bar 
has been developed to predict the SIP/RSI through-the-thickness interface 
stresses. The analysis assumes that the tile is rigid and that the SIP and 
filler bar behave as nonlinear spring foundations which are characterised by 
experimental stress-strain data. Interface stresses are calculated nua^rically 
from the tile equilibrium equations using a Newton iteration sches«. 

Parametric studies were conducted using this analysis for sianilated shock 
loading (combined load and moment) on square tiles. In addition, expected 
loads on an actual tile were also studied. The results of these studies indicate 
that when no mismatch is present, linear solutions tend to give nonconservative 
maximum stresses relative to the nonlinear solutions. When mismatch is present, 
the linear maximum stresses are only conservative at low load levels and highly 
nonconservative at high load levels. These studies also indicated that although 
the presence of filler bar reduces maximum tensile stresses in some cases, it 
does not appear that the nominal dimension filler bar, (1.27 cm or 1/2 inch wide) 
can reduce the nonlinear stress to a level at %ihich the linear solution becomes 
conservative relative to the nonlinear solution. Moreover, on an outboard 
trailing edge elevon tile at descf^' with * 2 ^ flapi (tile 139), the linear 
analysis is A9Z nonconservative relative to the nonlinear analysis with filler 
bar. 


INTRODUCTION 

Tests performed by Rockwell International and Langley Research Center 
indicate that the SIP possesses a highly nonlinear stress-strain relationship. 
It is the purpofe of this report to present an analysis to predict some of the 
possible consequences of the SIP nonlinearity on SIP/RSI interface through-the 
thickness stresses. Ihese stresses are considered since they appear to be 
critical to the integrity of the TPS. The application of this analysis to 
selected parametric studies is also presented and comparisons with a linear 
analysis are discussed. 
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GOVERNING EQUATIONS 
Assumptions 

Geometrical - In order to produce an efficient analysis applicable to a large 
number of tiles, SIP footprints and mismatch geometries, the following 
assumptions are made: 

(a) The tile behaves as a rigid body. 

(b) The tile rotations are small. 

(c) There exists an inherent mismatch between the tile and the substrate 
structure which is due to tile imperfections, tile warpage, substrate 
initial curvature or substrate deformation under load. 

Material - In order to account for the nonlinear SIP behavior observed experi- 
mentally, the following additional assumptions are made: 

(d) The SIP behaves as a nonlinear continuous spring - type foundation. 

(e) The stress strain relationship for the SIP may be accurately 
characterized as function of powers of the SIP strain. 


Equilibrium of Rigid Tile 


Consider the tile geometry shown in figures la, b where P is a positive 
vertical load acting at x = xp and y = yp, and and My are moments acting about 
the positive x and y axes respectively. Equilibrium of the tile results in 



where Aj^ and A^ are the portions of the tile 
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interface through-the-thickness stresses in the 
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As a consequence of assumptions (d) and (e). the constitutive relationship 
for the SIP may be expressed as 
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and are coefficients determined by experiment for 

>ar when j = 1 and for the SIP footprint when j » 2. Inasmuch a 


the filler bar 
the filler bar is not bonded to the tile. 

Furthermore, as a consequence of assumption (a) - (c) and figure 1 

W=Wo + '/‘^x -X«y 

where w is the local positive vertical tile displacement; w^ is the tile 
displacement at x = y == 0; and ta<y are rotations about the x and y 

axes respectively (positive as shown in figure la); and are the 

filler bar and SIP through-the-thickness strains respectively; t, and 
are the SIP and filler bar thicknesses respectively and, W is given by 


W = - Wj 

in which Wj is due to substrate deformation or mismatch and is due to 

tile imperfection or warpage. 

Substituting eqs. (4), (6), (7), and (8) into eqs, (1) - (3) yields 

^ (^0 - Ws + Wj:) c) - P = 0 

j Aj^ 

1 'X«<y-Ws+Wx'j dAj]-Mx-Pyp 
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Solution Procedure 

Equations (*>) - (11) may be solved for the three degrees of freedom, 

{ W., ). using a Newton iteration procedure. Thus the 

i * 1 / 

approximation to the solution vector is given by 
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where J is the symmetric Jacobian matrix 
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in which the residual vector is 
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and F, G, and H are the first terms in eqs. (9), (10) and (11) respectively. 
Also, 
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E, C"^. y) and are the filler bar and SIP tangent 

stiffnesses respectively which may be found by differentiating eq. (4) to give 



In order to provide the initial values of when 

r I in eq. (13), the following procedure comnonly employed in such non- 
linear problems is used. First, Pj and My in eqs. (9) - (11) are 

replaced by , and XMy respectively. The load parameter, \ , is 

then incremented from 0 to 1 in prescribed steps. At X = 0, the iterations 
of eq. (13) are initiated using W© = = 0 as trail values, and 

the iterations are allowed to converge. (These trial values are the exact 
solutions when no mismatch or substrate deformation is present; W - 0 ). 

At each succeeding value of X , the iteration procedure of eq. (13) is carried 
out using the converged solution at the preceeding value of X as the initial 
trail solution. 


. Uniqueness' of Solution 

It is shown in this section that for each set of prescribed stress-strain 
relationships governing the nonlinear SIP and filler bar behavior, only one 
solution of the equilibrium equations, eq. (1) - (3), exists independent of the 
order in which loads and/or substrate deformations are applied. To this end, 
the determinant of J is considered. 
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Expanding eq. (14) gives 

l^l = J^KyK<y +2 rMyKvy -WlJ Ky-my Ky-^^K-xy 

Since the axis system used for evaluating the eqs. (1) - (3) is arbitrary, a 

, -V A 

set of axes X ^ y may be sought about which 


my S Wy r Kxy =0 


The origin of this axis system V* ^ Yo relative to the x - y axes is 
given by 


-- iiijC yE;U,v)d4 

iii A * 


4 '* 

and its angular orientation, ^ , relative to the x *■ y axes is 


In the X - ^ coordinate system, 


1J|= ^KjK9 


The SIP and filler bar stress-strain relationships can be constructed such that 
£^^y^y^ is unique and positive, hence, 

■\ 7\>0 

Consider next the schematic load deflection curves in figure 7. If the 
dashed secondary solution curve exists, it must be a branch of the primary curve 
inasmuch as all solutions must yield zero load at zero deflection. However, since 
13| , there cannot be a branch (or bifurcation) point, and hence, only the 

primary solution exists. Consequently, once a set of stress-strain relationships 
for the SIP and filler bar are chosen, the solutions are'Ainique independent of 
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the order of load and/or substrate deformation application. However, to the 
extent that the selection of the stress-*strain relation depends o-i load history, 
the solutions will be load history dependent. 

Furthermore, since 1 3-) is always positive, Newton's iteration should 
converge for sufficiently small increments in X . Experience indicates 
that X increments of 0.1 are in most cases sufficient. 

DISCUSSION OF RESULTS 
Parametric Studies on Square Tiles 

A large number of tiles are of square planform with nominal dimensions of 
6" X 6" on a 5" x 5" SIP footprint. For the purpose of examining trends arising 
from parametric studies, it is useful to consider such tiles. (An actual tile 
of odd shape is considered at the end of this section). 

Effect of mismatch at different shock levels - Figure 2 presents the linear 
and nonlinear predicted maximum through-the-thickness SIP/RSI interface stress 
levels under a simulated transonic shock condition on a 3.81 cm (1.5 in) thick 
tile using Rockwell generated stress-strain data corresponding to 50 load cycles 
at +68.95 kPa (+10 psi^ as given in ref. 1. (The effect of other stress-strain 
data is considered later in this section). The tile is considered with and 
without a spherical mismatch of .048 cm (.019 in.) amplitude. Since the tile 
is considered rigid, the shock condition is statically equivalent to a combined 
load-moment condition or an eccentric load condition. 

The study indicates that under no external load,4p-0, the linear solution 
predicts higher stresses than the nonlinear analysis when mismatch is present. 

This is due to the very low stiffness (-about 27.58 kPa or 4 psi) at or near zero 

* 

load, associated with the nonlinear stress-strain data, while the linear solution 
uses a constant 172.4 kPa (25 psi) stiffness. However, as the shock amplitude is 
increased, the situation reverses and the nonlinear solution predicts considerably 
higher stresses than the linear solution; when A p 4,8 kPa (0.7 psi) the 

linear solution being up to 43% lower (nonconservative). The same tile without 
mismatch shows that the linear analysis is always nonconservative. Although not 
shown, these trends exhibited here are generally true for other loading conditions. 
Thus, due to the nonlinear SIP behavior, caution must be exercised in assessing 
the effect of mismatch, since as is shown, its affect varies with load level. 
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It is also of interest to consider the distribution of the SIP/RSI through*- 
the-thickness interface stresses. This is shown in figure 3. The figure shows 
that due to the low stiffness of the SIP at small strains, the majority of the 
tile carries relatively little stress, while the edges of the interface carry a 
highly concentrated stress. 

Effect of Filler Bar - Figure 4 shows the effect of different filler bar 
dimensions on the predicted nonlinear interface stresses for a square tile under 
combined load and moments. In this study, the SIP footprint dimensions are 
fixed Co be 12.7 x 12.7 cm^ (5x5 in^) and the surrounding filler bar dimension 
is increased from jt ^ ° 0 to 3.81 cm (1.5 in). Maximum interface stresses for 
different combinations of load and moment are *presented as ratios to the 
corresponding linear solution without filler bar, which is the value that the 
present Rockwell procedure would give for the same load and moment. Since the 
filler bar is only effective in compression, a moment must be applied to work 
the filler bar. For the study presented here that moment, M, was fixed at 16.9 
N*m (150 lb.' in.) and a direct tensile load, P, of variable magnitude was assumed 
to act simultaneously. 

The figure indicates that the nonlinear solution without filler bar = 0, 

is up to 70% higher than the RI solution and decreases with increasing in a 

nearly linear fashion for all values of P. As the load P is increased, the tile 
will be pulled further off of the filler bar which is only effective in compres- 
sion and, hence, the effect of the filler bar is decreased and the curves tend 
to become horizontal as the figure shows. However, for the nominal filler bar 
dimension used in the actual tiles, k. ^ ^ 1.27 cm (0.5 in), the presence of 
the filler bar is not sufficient to reduce the maximum stresses to the linear 
analysis value; the stresses still being 10 to 30% higher than the RI solution 
when J2. j = 1.27 cm (0.5 in). 

* 

Effect of Different Stress-Strain Data - As a consequence of the preceding 
discussion, it is clear that the linear solution may yield nonconservative 
results. However, as observed by others, the stress-strain behavior is strongly 
history dependent and exhibits considerable scatter from one specimen to another. 
The natural desire is to choose the stress-strain data giving the worst or most 
critical situation, however, this cannot be done without first examining the 
effect of the different stress-strain data. The purpose of this section is to 
pursue this examination on a square tile under shock loading. The stress-strain 
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data considered were generated by Rockwell and are as follows: 

(1) 2nd cycle data to +68.95 kPa (+10 psi); denoted as RI2 

(2) 50th cycle data to +68.95 kPa (+10 psi); denoted as RI50 

(3) 90Z Max data; denoted as RI90MAX 

(A) Average data; denoted as RIAV 

where (1) and (2) are given in reference 1, and (3) and (A) are given in 
reference 2. For all data, the stress-strain curves are linearly extrapolated 
above the strain levels for which the data are valid. 

A comparison of the results shown in figure 5 reveals a relatively wide 
variation when different stress-strain data is used in the presence of mismatch, 
however, the trends discussed previously for the RI50 data remain unchanged. It 
may be argued that the RI50 results are too conservative since their use implies 
that the tile has seen 50 cycles at 68.95 kPa (^10 psi) which is quite high. In 
fact, high stress cycles tend to increase the deadband which, in turn, produces 
higher stress concentrations at tile corners while cycles at lower stress levels 
do little to extend the deadband. Therefore, the R12 curve is probably more 
reasonable. 


Outboard Elevon Trailing Edge Upper Surface Tile 

In addition to the parametric studies described previously, the analysis has 
also been appUed to the outboard elevon trailing edge upper surface tile, #139, 
as shown in tigure 6 for descent at +20^ flap using RI2 data extended to +103. A 
kPa (+15 psi) from curve 6M2A of ref. 3. The tigure illustrates the highly non- 
linear interface stress distribution predicted on such an odd shaped tile under 
the expected load and moments. For the purpose of comparing the linear and non- 
linear solutions, a linear solution without filler bar was used since this 
corresponds to the Rockwell analysis procedure. 

From the figure it is observed that the maximum nonlinear stress when filler 
bar is absent is about 53% larger than the maximum linear stress. When fill.jr 
bar is added, this percentage decreases slightly to A9%. The addition of the 
filler bar helps to distribute the compressive stresses so that these are highly 
reduced, but has little effect in reducing the tensile stresses. 
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CONCLUSIONS 

A nonlinear analysis has been developed to predict the SIP/RSI through-the- 
thickness interface stresses including the effects of mismatch and filler bar. 
The analysis assumes that the tile is rigid and the SIP and filler bar behave 
as a nonlinear spring foundation which is characterized by experimental stress- 
strain data. Interface stresses are calculated numerically from the equilibrium 
of the tile using a Newton iteration scheme. • 

As a result of parametric studies for a simulated shock loading on square 
tiles, the following conclusions are derived: 

1. When no mismatch is present, linear solutions for maximum interface 
stresses tend to be nonconservative relative to nonlinear solutions. 

2. When mismatch is present, linear solution relative to nonlinear 
solutions tend to be only conservative at low load levels and become 
highly nonc( aservative at higher levels. 

3. Use of filler bar reduces no. linear stresses, but for nominal 
dimension filler bar up to 1.27 cm (V) wide, reductions are 
generally small, thus linear analysis without filler bar as used by 
Rockwell still underestimates the nonlinear analysis with filler bar. 

A. Conclusions 1, 2, and 3 are generally independent of the choice 

of different Rockwell stress-strain data; however, the magnitude of 
the difference between linear and nonlinear solutions can vary 
significantly when large mismatches are present. 

Tlie analysis has also been applied to the outboard trailing edge elevon tile, 
(tile 139), for descent at •t-20‘^ flap. Results for this tile indicate that the 
linear analysis underestimates the maximum interface stress by 49% even when 
filler bar is present. 
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Figure 1. TPS Geometry: 

(a) SIP footprint and filler bar planform 

(b) RSI tile profile 
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Figure 2. Variation of maximum interface stress with simulated 
shock magnitude for a square tile with and without 
spherical mismatch. 
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Stress at y'*0 for 20.7 kPa (3 PSI) Simulated Transonic 
Shock for RI50 and Linear (E=17.2 kPa=25 PSI) SIP on 
Tile with Spherical Mismatch. 
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Figure 5. Effect of Different Stress-Strain Data on Maximum 
Through- the-Thickness SIP/RSI Interface Stress for 
Shock on Square Tile With Spherical Mismatch 
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Figure 7. Schematic of possible load-deflection paths 


